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Abstract 

The momentum transfer dependence of the total cross section for elastic scattering of 
cold dark matter candidates, i.e. lightest supersymmetric particle (LSP), with nuclei is 
examined. The presented calculations of the event rates refer to a number of representative 
nuclear targets throughout the periodic table and have been obtained in a relatively wide 
phenomenologically allowed SUSY parameter space. For the coherent cross sections it is 
shown that, since the momentum transfer can be quite big for large mass of the LSP and 
heavy nuclei even though the energy transfer is small (< lOOiLeU), the total cross section 
can in such instances be reduced by a factor of about five. For the spin induced cross 
section of odd-A nuclear targets, as is the case of studied in this work, we found 

that the reduction is less pronounced, since the high multipoles tend to enhance the cross 
section as the momentum transfer increases (for LSP mass < 200 GeV) and partially cancell 
the momentum retardation. The effect of the Earth’s revolution around the sun on these 
event rates is also studied by folding with a Maxwellian LSP-velocity distribution which is 
consistent with its density in the halos. We thus found that the convoluted event rates do 
not appreciably change compared to those obtained with an average velocity. The event 
rates increase with A and, in the SUSY parameter space considered, they can reach values 
up to 140 y~^Kg~^ for Pb. The modulation effect, however, was found to be small (less than 
±5%). 
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I. INTRODUCTION 


There is ample evidence that about 90% of the matter of the universe is dark i-i- There are, 
also, numerus arguments indicating that our galaxy is immersed in a dark halo which outweights 
the luminous component by a factor of about ten. Furthermore, the large scale structure of the 
universe, may be accommodated supposing two kinds of dark matter |^. One kind is composed 
of particles which were relativistic at the time of the structure formation. This is called Hot Dark 
Matter (HDM). The other kind is composed of particles which were non-relativistic at the time of 
structure formation. These constitute the Cold Dark Matter (CDM) component of the universe. 
In any case the CDM component of the universe is at least 60%. Obviously, the galactic halo is 
composed of dark matter, since HDM particles will be moving too fast to be trapped in galaxy. 
There are two candidates for CDM. The hrst is MACHO’s (Massive Compact Halo Objects, 
i.e. white dwarfs, Juppiter-like objects etc.) and the other is more exotic, i.e. WIMP’s (Weak 
Interacting Massive Particles). Recent preliminary experiments [Q suggest that about half of the 
mass of the halo is made of WIMP’s. The most appealing possibility linked with supersymmetry 
(SUSY) of a WIMP candidate is the LSP, i.e. the Lightest Supersymmetric Particle ||]-[0 (see 
Ref. 0 for a recent review). 

In recent years, the phenomenological implications of Supersymmetry are being taken very 
seriously More or less, accurate predictions at low energies are now feasible in terms of 

few input parameters in the context of SUSY models P]-[|l. Such predictions do not appear to 
depend on arbitrary choices of the relevant parameters or untested assumptions. In any case the 
SUSY parameter space is somewhat restricted ^-[0. 

In such theories derived from Supergravity the LSP is expected to be a neutral Majorana 
fermion with mass in the 10 — 500GeU/c^ region travelling with non-relativistic velocities (< 
(3 >~ 10“^), i.e. with energies in the KeV region. In practice, however, one expects a velocity 
distribution which is supposed to be Maxwellian (see sect. IV). In the absence of R-parity violating 
interactions this particle is absolutely stable. But, even in the presence of R-parity violation, it 
may live long enough to be a CDM candidate. 

The detection of the LSP, which is going to be denoted by Xi, is extremely difficult, since 
this particle interacts with matter extremely weakly. One possibility is the detection of secondary 
high energy neutrinos which are produced by pair annihilation in the sun where this particle is 
trapped. Such high energy neutrinos can be detected via neutrino telescopes. 

The other possibility, to be examined in the present work, is the detection of the energy of the 
recoiling nucleus {A, Z) in the reaction 

Xl + (^) Z) ^ Xl + zy ( 1 ) 

This energy can be converted into phonon energy and detected by a temperature rise in cryostatic 
detector with sufficiently high Debye temperature P, |^, |^. The detector should be large enough 
to allow a sufficient number of counts but not too large to permit anticoincidence shielding to 
reduce background. A compromise of about IKg is achieved. 

Another possibility is the use of superconducting granules suspended in a magnetic field. The 
heat produced will destroy the superconductor and one can detect the resulting magnetic flux. 
Again a target of about IKg is favored. 

There are many targets which can be employed. The most popular ones contain the nuclei 
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iJ/e, fjNa, »Si, ‘«Ca, ^^Ge, Jf/, Jf A'e and ^Pb. 

In order to be able to calculate the event rate for the process (1) the following ingredients are 
necessary. 

1) One must be able to construct the effective Lagrangian at the elementary particle level in 
the framework of Supersymmetry |^-|^. We will follow the procedure adopted in Ref. [ 0 - For 
the readers convenience we will provide the important elements in sect. II. 


2 ) One must make the transition from the quark to the nucleon level This is not 


straightforward for the scalar couplings, which dominate the coherent part of the cross section, 
and the isoscalar axial current which is important for the incoherent cross section for odd targets. 


3) One must properly treat the nucleus. Admittedly, the uncertainties here are smaller than 
those of even the most restricted SUSY parameter space. One, however, would like to put as 
accurate nuclear physics input as possible in order to constrain the SUSY parameters as much as 
possible when the data become available. 


For the coherent production, if one ignores the momentum transfer dependence, the proce¬ 
dure is straightforward. The spin matrix element, however, is another story. For its evaluation 
practically every known nuclear model has been employed. 


At hrst, the Independent Single Particle Shell Model (ISPSM) had been employed 0 0, |5|, 
Subsequent calculations using the Odd Group Method (OGM) and the Extended Odd Group 
Method (EOGM), utilizing magnetic moments and mirror /d-decays, by Engel and Vogel 0], 
showed that the ISPSM was inadequate (see also Ref. [|^). Eventually, however, by performing 
shell model calculations 0 , this model was also found lacking (see Ref. p9[). lachello, Krauss and 
Maino m employed the Interacting Boson Fermion Model (IBFM) and Nikolaev and Klapdor- 
Kleingrothaus HI. 0 the finite fermion theory in order to reliably evaluate the spin matrix 
elements. 


One additional complication arise from the fact that the LSP appears to be quite massive, per¬ 
haps heavier than lOOGeU. For such heavy LSP and sufficiently heavy nuclei, the dependence of 
the nuclear matrix elements on the momentum transfer cannot be ignored even if the LSP has en¬ 
ergies as low as lOOiFeU) |0, |^. This affects both the coherent and the spin matrix elements. 

For the coherent mode the essential new features can be absorbed in the nuclear form factor. The 
evaluation of the spin matrix elements is quite a bit more complicated. Quite a number of high 
multipoles can now contribute, some of them getting contributions from components of the wave 
function which do not contribute in the static limit (i.e. at q=0, see sect. III). Thus, in general, 
sophisticated Shell Model calculations are needed to account both for the observed retardation of 
the static spin matrix element and its correct q-dependence. for the experimentally interesting 
nuclear systems ^Si and l^Ge very elaborate calculations have been performed by Ressell et al. 
In the case of 32^6 a further improved calculation by Dimitrov, Engel and Pittel has recently 
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been performed by suitably mixing variationally determined triaxial Slatter determinants. 
Indeed, for this complex nucleus many multipoles contribute and the needed calculations involve 
techniques which are extremely sophisticated. 


From the above discussion the necessity for more detailed calculations especially for the spin 
component of the cross sections for heavy nuclei is motivated. The aim of the present paper is 
to calculate LSP-nucleus scattering cross section using some representative input in the restricted 
SUSY parameter space as outlined above. The coherent matrix elements are computed throughout 
the periodic table (sect. III). The needed form factors were obtained using the method of Ref. 
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34|, which are in good agreement with experiment. For the spin matrix elements we have chosen 


^Ph as target. This target, in addition to its experimental qualifications, has the advantage of a 
rather simple nuclear structure (see sect. III). Thus, only two multipoles can contribute. 


Finally, the counting rates are folded with a reasonable LSP-velocity distribution 0 (see sect. 
IV) in order to estimate the modulation due to the Earth’s velocity. Convoluted rates are obtained 
(sect. IV) in the framework of models discussed in the introduction and sect. II. 


II. BRIEF DESCRIPTION OF THE OPERATORS 


It has recently been shown that process (|l|) can be described by a four fermion interaction 


of the type |]n 


^e// - - 




[^AXidVxi + Jxixi] 


( 2 ) 


where 


J\ = N'^x[fv + fv^3 + if A + /A 7 - 3)75 ]N 


( 3 ) 


and 


J = Nifl + fsr3)N 


( 4 ) 


We have neglected the uninteresting pseudoscalar and tensor currents. Note that, due to the 
Majorana nature of the LSP, = 0 (identically). The vector and axial vector form factors 

can arise out of Z-exchange and s-quark exchange |]^-|^ (s-quarks are the SUSY partners of 
quarks with spin zero). They have uncertainties in them (for three choices in the allowed parameter 
space of Ref. [§ see Ref. |j^). In our choice of the parameters the LSP is mostly a gaugino. 
Thus, the Z- contribution is small. It may become dominant in models in which the LSP happens 
to be primarily a Higgsino. Such a possibility will be examined elsewhere. The transition from 
the quark to the nucleon level is pretty straightforward in this case. This is in general the case of 
vector current contribution. We will see later that, due to the Majorana nature of the LSP, the 
contribution of the vector current, which can lead to a coherent effect of all nucleons, is suppressed 
[10-16]. The vector current is effectively multiplied by a factor oi f3 = v/c, v is the velocity of 
LSP (see Table I). Thus, the axial current, especially in the case of light and medium mass nuclei, 
cannot be ignored. 


For the isovector axial current one is pretty confident about how to go from the quark to the 
nucleon level. We know from ordinary weak decays that the coupling merely gets renormalized 
from qa = 1 to qa = 1.24. For the isoscalar axial current the situation is not completely clear. 
The naive quark model (NQM) would give a renormalization parameter of unity (the same as the 
isovector vector current). This point of view has, however, changed in recent years due to the 
so-called spin crisis, i.e. the fact that in the EMC data |T^ it appears that only a small fraction 
of the proton spin arises from the quarks. Thus, one may have to renormalize by = 0.28, 
for u and d quarks, and = —0.16 for the strange quarks [|^], i.e. a total factor of 0.12. These 


two possibilities, labeled as NQM and EMC, are listed in Table 1. One cannot completely rule out 
the possibility that the actual value maybe anywhere in the above mentioned region . 
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The scalar form factors arise out of the Higgs exchange or via s-quark exchange when there 
is mixing between s-quarks ql and qR []^-|^ (the partners of the left-handed and right-handed 
quarks). They have two types of uncertainties in them. One, which is the most important, at the 
quark level due to the uncertainties in the Higgs sector. The other in going from the quark to the 
nucleon level O, M. Such couplings are proportional to the quark masses, and hence sensitive 


to the small admixtures of qq (q other than u and d) present in the nucleon. Again values of fg 


and fg in the allowed SUSY parameter space are considered 0- 

The actual values of the parameters and used here, arising mainly from Higgs exchange, 
were obtained by considering 1-loop corrections in the Higgs sector. As a result, the lightest Higgs 
mass is now a bit higher, i.e. more massive than the value of the Z-boson pT| , The thus 
obtained values of the parameters /_§ and fg are smaller than those of Ref. [|^ (see Table I). The 
next source of ambiquities involves the step of going from the quark to the nucleon level for the 
scalar and isoscalar couplings. Here we adopt the procedure described in Ref. [|^ 
the analysis of Ref. O, |2f 


as a result of 


III. TOTAL CROSS SECTION 


The invariant amplitude in the case of non-relativistic LSP takes the form 


mf 

~ /?Vor+iJp+i^r 


(5) 

where mi is the LSP mass, Jo and J indicate the matrix elements of the time and space 
components of the current Jx of Eq. (§), respectively, and J represents the matrix element of 
the scalar current J of Eq. (H). Notice that Jq P is multiplied by (5"^ (the suppression due to the 
Majorana nature of LSP mentioned above). It is straightforward to show that 

\j»? - A^\F(e)\‘ [fv 

(6) 

= /l=|F(q=)r (/S 

(7) 

|jp = < J.ll [/On„(q) + /itJi(q)] ||J. > 

(8) 


with E(q^) the nuclear form factor and 


4 4 

no(q) = Z n,(q) = aU)T,U)e-‘‘’‘- (9) 

i=i i=i 

where cr(j), are the spin, third component of isospin (rsjp >= \p >) and coordinate of 

the j-th nucleon and q is the momentum transferred to the nucleus. 
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The differential cross section in the laboratory frame takes the form [17 


da aQ nil 2 

dVt 71 rrip (1 + rfY 




(l + r()= 


( 10 ) 


where rrip is the proton mass, r] = mi/nipA, ^ = Pi ■ 0 (forward scattering) and 


(Jo = —{GfitIp) ~ 0.77 X 10 

ZTT 


( 11 ) 


The momentnm transfer q is given by 


|q| = qo = P 


2miC 
I+7] 


( 12 ) 


Some values of qq (forward momentum transfer) for some characteristic values of mi and 
representative nuclear systems (light, medium and heavy) are given in Table II. It is clear from 


Eq. (h^ that the momentum transfer can be sizable for large mi and heavy nuclei {rf small). 
The total cross section can be cast in the form 


a = 


+ 

+ 


o'o( 


mi 2 


m„ 


{i+vy 


{71^ 


A-2Z 

A 




(fs fs 


lA 


2\ + 1 Q ^A 

J-'oWoJ —T7T7— ^\JV ~ JV 


2Z^ 




^ 2 J-UV.U2 2(1+77)2^^" A 

(fA^o(0))^Ioo(qo) + 2/A/Allo(0)fli(0)/oi(go) + (/Alli(0))^-^ii('?o) } 


(13) 


The quantities entering eq. (0) are explained in the subsections A and B below. 

A. The coherent matrix element 


The terms in the square bracket of Eq. ([TB|), describe the coherent cross section the momentum 
dependence of which is involved in the integrals Ip{qQ), where p =0 for the isoscalar and p = 1 for 
the isovector component. These integrals are given by 


/,(<(o)=2(p+i) !'e^^'’\F(qt,e)?di, p=o,i 
0 


(14) 


The integrals Ip are normalized so as lp(0) = 1 (following the normalization E(0) = 1) and 
they can be calculated by using Ref. [^] where we have shown that, the nuclear form factor F{q^) 
can be adequately described within the harmonic oscillator model as 


7 N 

F{q^) = [ jHqb, Z) + -$(g6 , N) ] 

(N = A — Z) where b ~ lAA^A and <h is a polynomial of the form 


( 15 ) 
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a = Z,N. 


(16) 


Nmax (ci) 

$(g6,a)= ^ d^^\qbf^, 

A=0 

Nmax{Z) and Nmax{N) depend on the major harmonic oscillator shell occupied by protons and 
neutrons respectively. The integrals /p(q'o) can be written as 

Ip{qo) = Ip{u) = (1 + f \F{2x/b^)\‘^ dx, (17) 

0 

where 

u = qlby2 (18) 

With the use of Eqs. (H), (0) we obtain 

/,(«) = T{ + 2NZli;l(u) + /VV<fk«)} (19) 

where 




Nmax{a) Nmaxi/S) q{oi) ^(/3) ^X + U+p _j_ ^ _j_ p^j 


= (i+rt i: i: 


a (3 


u 


i+p 


X+V+p K 

^ k\ J 


( 20 ) 


(a,/? = N^Z). The coefficients are given in Ref. [^], for light and medium nuclei, and in 
Ref. 1^ for heavy nuclei. 


The integrals Jq for three typical nuclei and "^Pb ) are presented as a function of 

rrii in Fig. 1. The values of the harmonic oscillator parameter b used are: b=1.849 fm for Ca, 
b=2.039 fm for Ge, and b=2.434 fm for Pb. We see that, for light nuclei the modihcation of the 
cross section by the inclusion of the form factor is small. For heavy nuclei and massive mi the 
form factor has a dramatic effect on the cross section and may decrease it by a factor of about 
Eve. The integral R is even more suppressed but it is less important (see Ref. USSj). 


B. The spin matrix element 


The other terms in Eq. (pRf), describe the spin dependence of the cross section in terms of 
the Ipp/, with p, p' = 0,1. The latter integrals result by following the standard procedure of the 
multipole expansion of the in Eq. (P) and are dehned by 


^pp' i%) 



SJ,(0) 0/(0) 


p, p =0,1 


where we have made the identihcation 


( 21 ) 


Slj,”'" = = (2Ji + 1) ) < Jf\\J23o(.1a^rj)uii,U)a(j)\\J, >, /3 = 0,1 (22) 




with uJoU) = 1 = 'TsU)- general. 
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^(^Aa) = (2 j. + 1) l^/^ < Jj\\Y^\^^{rj)®a{j)Yix{qoirj)u;p{j)\\Ji> p = 0,1 (23) 


i=i 


With the above expressions we have managed to separate the elementary parameters f\ and f\ 
from the nnclear parameters. 


We warn the reader that, the integrals of Eq. (|^) are normalized to nnity as g —>■ 0, 


i.e. 


Ipp'{qo = 0) = 1. This normalization is different than that fonnd in the previons literatnre. In the 
above limit (go = 0) the spin matrix element takes the simple expression 


J 


2 


/“SJ„(0) + /‘SJi(O) 


The spin matrix elements, nnfortnnately, depend in general rather sensitively on the details of 
the nnclear structnre, which is inclnded in the integrals /pp/(go) (see Eq. ( 0 ). As we mentioned 
in the introduction, the hrst attempt for quantitative description was based on the Odd Group 
Method 1^. Subsequent shell model calculations demonstrated that the OGM was not adequate 
and showed that more elaborate calculations were needed |^, Furthermore, since the matrix 
elements at g = 0 are often quenched, the momentum dependence of the matrix elements was 
more important than it was naively expected. As a matter of fact, one has to include a lot of 
conhgurations to accommodate all multipoles, which in complex nuclei like and result 

in very large Hilbert spaces. It will be therefore, a very hard task to substantially improve the 
calculations of Ref. 


32-33 


Among the targets which have been considered for LSP detection, stands out as an 

important candidate. The spin matrix element of this nucleus has not been evaluated, since one 
expects the relative importance of the spin versus the coherent mode to be more important on light 
nuclei. But, as we have mentioned, the spin matrix element in the light systems is quenched. On 
the other hand, the spin matrix element of especially the isoscalar one, does not suffer from 

unusually large quenching, as is known from the study of the magnetic moment. Thus, we view it 
as a good theoretical laboratory since: i) It is believed to have simple structure, one 2pl/2 neutron 
hole outside the doubly magic nucleus ii) Because of its low angular momentum, only two 

multipoles A = 0 and A = 2 with a J-rank of k = 1 can contribute even at large momentum 
transfers. One can thus view the information obtained from this simple nucleus as complementary 
to that of which has very complex nuclear structure. 


To a good approximation the ground state of the Pb nucleus can be described as a 2pi/2 


neutron hole in the l^^Pb closed shell. Then for A = 0 one hnds 


Go(q) = -(l/v^)P 2 p(q'), Gi(q) = {l/V3)F2p{q^) (24) 

and 


Joo = /oi = /ii = 2 e [F 2 p{q^)? di (25) 

Even though the probability of Ending a pure 2pi/2 neutron hole in the 2 ground state of ^Pb 
is greater than 95%, the ground state magnetic moment is quenched due to the I’*' p-h excitation 
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involving the spin orbit partners. Hence, we expect a similar suppression of the isovector spin 
matrix elements. Thus we write 


1 ( 1 / 2 )- ><„ 


— Col(2^1/2) ^ > 


1 + Ci| [0iii/2(/^)(0il3/2) > 


+ C 2 I[0h9/2(p)(0hii/2) > +••• 


(26) 


Due to angular momentum and parity selection rules, we have k = 1 and A = 0, 2. Retaining 
terms which are at most linear in the coefficients Ci, C 2 we obtain 

i) A = 0 


Do(q) = Cl {F 2 p(g 2 )/v^ - 8 [(7/13)iC^iFo*(g') + ( 5 /ll)'C^ 2 Co,(g 2 )] } (27) 

Di(q) = - Cl {F2p(g')/V3 - 8 [(7/13)'CCiFo*(g') - {h/llf'^C 2 Foh{q^)] } (28) 

where 

^max 

Fniiq^) = ( 29 ) 

/i=0 

and the coefficients Cq, C\ and C 2 were obtained by 
, ^ in a model space of 2h-lp conhgurations. They 

are given by 
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The coefficients 7 )”'^ are given in Ref. 
diagonalizing the Kuo-Brown G-matrix 


Co = 0.973350, Ci = 0.005295, C 2 = -0.006984 

We also hnd 

Do(0) = —(1/a/3) (0.95659), {small quenching) (30) 

f^i(O) = (l/v^)(0.83296), {sizable quenching) (31) 

The amount of retardation of the total matrix element depends on the values of and f\. 
ii) A = 2. 

In this case, in addition to the leading ( 2 ^ 1 / 2 )”^ conhguration, the hrst leading correction to 
the nuclear matrix element is linear in the mixing coefficients Cj^jj appearing in the expression: 

l(^)”^ >= Co{|(2pi/2)"^(n) > +^Cj,j2|(2pi/2)"^(n); {ji^j2)Ji2 = 1; ^ > } (32) 

hh 


i.e. 












(33) 


^ V2 J, + 1 


{2 < nih\j2{qo^r)\n2l2 > 


J1J2 


+{-ir < {2p,/2)-m\i^pi/2)-^ >} 


where G{ji,j 2 ,p) are isospin dependent geometrical factors which can be evaluated by standard 
techniques. The radial integrals < ni/i|j 2 (<?o^^)|^ 2^2 > can be cast in the form of Eq. (|^) , but 
they depend on two single particle quantum numbers (see appendix). The relevant coefficients 
used in the present work are given in Ref. . 

Notice, however, that unlike the A = 0 case, many amplitudes can contribute if the quadrupole 
modes coupled to the single hole wavefunction are admixed in the ground state of the system. 


In the simple model of Ref. [0, in addition to the Ci, C 2 encountered above, one needs 


the amplitudes of the two additional configurations, 0 ji 3 / 20 jii/ 2 (n) and ^iii/ 2 ^g 9 / 2 {p)i which are 
Ca = 0.000239 and C 4 = —0.000642. Obviously, this is a simplification, since one should consider 
the spin Giant Quadrupole Resonance (GQR), which may have a small admixture in the ground 
state of the nucleus but a very large transition matrix element. Such a detailed calculation 
including all 2tiuj excitations is in progress and it will be reported elsewhere. 


Using Eqs. (1^) , (^), (^ 8 [) and (|33|) , we can evaluate the integrals /qo, /qi and Jqo for 
The results for In are presented in Fig. 2 (the other two are practically indistinguishable, see Ref. 
|]35|). In Fig. 2(a) In is plotted as a function of the LSP mass while in Fig. 2(b) it is plotted 
(together with Jq) as a function of the parameter u given in Eq. (|^. We see that for a heavy 
nucleus and high LSP mass the momentum transfer dependence of the spin monopole (A = 0) 
matrix elements is quite large. It is, however, to a large extent neutralized by the spin quadrupole 
(A = 2). So the overall effect is not dramatic for LSP mass less than lOOGeU. As a matter of 
fact, from Fig. 2(b) we see that, the retardation of the spin matrix element is quite a bit less than 
that of the coherent mode for almost all values of u. 


IV. CONVOLUTION OF THE CROSS SECTION WITH 
THE VELOCITY DISTRIBUTION 


The cross sections which would be given from an LSP-detector participating in the revolution 
of the earth around the sun would appear retarded. In this section we are going to study this 
effect by using the method of folding. To this aim let us assume that the LSP is moving with 
velocity Vz with respect to the detecting apparatus. Then the detection rate for a target with 
mass m is given by 


dN 

dt 


p( 0 ) m 
rrti Arrip 


Vz\a{v) 


(34) 


where p(0) = 0.3GeU/cm^ is the LSP density in our vicinity. This density has to be consistent 
with the LSP velocity distribution. Such a consistent choise can be a Maxwell distribution 


fiv') = (i/ift>o) “e I'''''™)’ 


(35) 


provided that 


H 
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Vo — \j (2/3) < v'^ > — 220 Km /3 


(36) 


For our purposes it is convenient to express the above distribution in the laboratory frame, i.e. 

/(v.VE) = (V)?«(,)-=e-<'’+''*>"'”« (37) 

where v^; is the velocity of the earth with respect to the center of the distribution. Choosing a 
coordinate system in which X 2 is the axis of the galaxy, X 3 is along the sun’s direction of motion 
(vq) and xi = X2 X X3, we find that the position of the axis of the ecliptic is determined by the 
angle 7 ~ 29.80 (galactic latitude) and the azymouthal angle uj = 186.3° measured on the galactic 
plane from the X 3 axis [^ . 

Thus, the axis of the ecliptic lies very close to the X 2 X 5 plane and the velocity of the earth is 


ve = Vq + vi = Vq + Ui (sma xi — cosa cosy X 2 + cosa smy X 3 ) (38) 


and 


cos a 


Vo ■ Vi = Uo^’l 


\/l + cot'^'j COs'^OJ 


vqVi sin y cos a 


(39) 


where vq is the velocity of the sun around the center of the galaxy, v\ is the speed of the earth’s 
revolution around the sun, a is the phase of the earth orbital motion, a = 2n{t — ti)/TE, where 
ti is around second of June and Te = lyear. 

The mean value of the event rate of Eq. (^^ , is defined by 



p( 0 ) m 
mi Amp 


/(v, ve) I I a(|v|)(i°v 


Then we can write the counting rate as 


(40) 


p(0) m I -T— 


where 


(41) 


< S >= 


V< > 

Thus, taking the polar axis in the direction ve, we get 


/(v, VE)cr(|v|)ci°v 


< S >= 




r v^dv j' |e|dee-<”‘+''n2™vS)/”2c(ti) 


(42) 


(43) 


or 
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< s >= 


2 




vdv Fq 


‘2'VVE-. _(y2: 2),2 . , 

n 


(44) 


with 


^o(x) = Xsinhx - coshx + 1 
One can also write Eq. (^41) as followes 


,0, 1 /■OO q, 

< S >= (-) M — fi{v)a{v)dv 
^ 3 ^ Jo no 


(45) 


(46) 


with 




1 


In the case in which the hrst term in Eq. (^) becomes dominant, we get 


(47) 


hiy) 



{v - Ve^I 
.,,2 


— exp 


{v + VEfl 'I 


(48) 


in agreement with Eq. (8.15) of Ref. |]^. In Eq. (|4^ the nnclear parameters are implicit in 
the cross section a(v) given from Eq. (|T^) . The nnclear physics dependence of < S > could be 
disentangled by taking note of the extra velocity dependence of the coherent vector contribution 
in a(v) and introducing the parameters 


6 = —— = 0.27, V’ = —) (49) 

where the quantity uq is the one entering the nuclear form factors of Eq. (O) for V = Vo, which 
in this case is given by 


Wo — X 


1 /2/domic^ b 


Afterwards, we can write Eq. 


as 


2 I (1 + r;) he 


Po = - 
c 


(50) 


< E > = 


mi \ 2 (To 


(-) 

\rn. / 


(1 + ?7)2 




,A-2Z 

Jv 




(51) 




A ' 2(l + i;)2 

'(/.4^o(0)) Joo + 2/)[/^f2o(0)fli(0)Joi + (^/^f2i(0)^ Juj 


lA 2Z 2 j 

—) Jo 
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If we assume that Jqo = >^01 = Jiii as seems to be the case for the spin dependent part 

r- 1 2 

of Eq. (^7|) is reduced to the familiar expression /°f2o(0) + /^f2i(0) Jn, where the quantity in 
the bracket represents the spin matrix element at g = 0. 


The parameters Jq, Jp, Jpa describe the scalar, vector and spin part of the counting rate, 
respectively, and they are given by 


2 e~^ 

Jo(A,no) = ^=-— / Fo{2\-il))Io{uQ'ip‘^)d'il) (52) 

VbTT Jo 

Jp{\uo) = —Fo(2Ai/’)/p(Mo'0^)# (53) 

V67r Jo 

2 r°° ;2 

JpaiKuo) = -j=—— / FQ{2\'il))Ip„{uoi)‘^)dil) (54) 

VbTT A^ Jo 

1) TP V o1 1/2 

A = —= 1 F dcosasin'-^ + (5/ 2 Y (55) 

Vo L J 


The parameters Jp, Ip^ have been discussed in the previous section. The above integrals are 
functions of A and uo- The latter depends on vq, the nuclear parameters and the LSP mass. These 
integrals can only be done numerically. Since, however, A is close to unity, we can expand in 
powers of 5 and make explicit the dependence of these integrals on the earth’s motion. Thus, 


~ ^ r ~ ~ 

Jo(A,mo) = — Bi Kq°\uo) + d sinj cosa Kq^\uo) 

vGtt 


(56) 


2 

Jp(\uo) ^ -^B2 + i^3mjcosaAF(uo) 


(57) 


2 

J ^ 0 ) [A'iS’(“o) +i5sM!7cosaii'™(«o) 


(58) 


The integrals Kq, and are normalized so that they become unity at mq = 0 (negligible 
momentum transfer). We hnd 





Fo{2p)dp 


/ + 2z/ 1.860 

e 


(59) 


2 9 2 S 

B 2 = — Foi2p)dp = -( - + 7z/) ^ 4.220 

de Jo 3 e 

with 


(60) 


z/ 


e ^^dt 


0.747 


(61) 
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Furthermore, 



Ko = ^ 

° eBi 

roo ^ 

/ Fi{2'ilj)Io{uoiJ 

Jo 

Z = 0,1 

(62) 


K‘= 2 

" 36^2 

poo 2 

/ Fi{2ij)Ip{uoij^)dij, 

Jo 

Z = 0,1 

(63) 


e5i 

POO 2 

/ '06“’^ Fi{2'ijj)Ipkuo'ilJ )d'ilj, 

Jo 

Z = 0,1 

(64) 

with Fq{x) given in Eq. (^) and 





Fdx) ■ 

r 

= 2 [(—+ l)coshx - ysmhx 

- l' 

(65) 

The counting rate 

can thus be cast in the form 





/dN \ /dN\ 

(if) = (if)X + 


(66) 


where (^)o is the rate obtained from the I = 0 multipole and h the amplitude of the oscillation, 
i.e. the ratio of the component of the multipole Z = 1 to that of the multipole I = 0. Below (see 
also Tables III and IV) we compute separately the amplitude of oscillation for the scalar, vector 
and spin parts of the event rate i.e. the quantity h = 5siw^ {uq)/{uq) . Note the presence of 
the geometric factor siny = 1/2, which reduces the modulation effect. 

In order to get some idea of the dependence of the counting rate on the earth’s motion, we 
will evaluate the above expressions at Uq = 0. We get 


^ 0 ° = K% = I 


kl = Kl^ = — - 

° l/e + 2i/ 


0.402 


^ 3/(2.)+(ll/2)./2 ^ 
^ 3/e + 7z/ 


Thus, for siny ~ 0.5 


2 


Jq ^ 1.860(1 + 0.054cosa) = 0.857(1 + O.Obicosa) 


2 


Jp = ^=4.220(1 + O.OOOcosa) = 1.944(1 + 0.099cosa) 


(67) 

( 68 ) 

(69) 

(70) 

(71) 
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We see that, the modulation of the detection rate due to the earth’s motion is quite small 
{h ~ 0.05). The corresponding amplitude of oscillation in the coherent vector contribution, Eq. 

is a bit bigger {h ^ 0.10). However, this contribution is suppressed due to the Majorana 
nature of LSP (through the factor /3^). The modulation due to the Earth’s rotation is expected 
to be even smaller. 

The exact integrals, for the 1=0 and 1=1, are shown in Figs. 3(a)-3(c). The most important 
of these integrals, those of Eq. (|6^ associated with the scalar interaction, are shown in Fig. 3(a). 
In Fig. 3(b) we present the integrals of Eq. ([H3| ) for p = 0 associated with the vector interaction 
(the integral for p = 1 is analogous but it is less important). Finally in Fig. 3(c) the integrals of 
Eq. ([6^) for p = 1 and <7 = 1 are shown. The others are practically indistinquishable from these 
and are not shown (see Ref. |]35|). 

Before closing this section we should mention that, the folding procedure can also be applied 
in the differential rate in order to obtain the corresponding convoluted expression for da/dQ, i.e. 
before doing the angular integration in Eq. (pTf ) and obtain the total cross section Eq. (pRj) . The 
resulting expressions are, however, a bit more complicated and they will not be given here. 

V. RESULTS AND DISCUSSION 


The main goal of the present work was the calculation of the cross sections for the scattering of 
LSP, a Cold Dark Matter Candidate, with nuclei. The coherent scattering was evaluated for three 
typical nuclei, Ca, Ge and Pb, and the spin matrix element for The momentum dependence 

of the matrix elements was taken into account. Special attention has been paid to evaluating the 
modulation of the event rates due to the earth’s motion. Our results are summarized as follows. 

A. Coherent scattering 


As we have seen in sect. Ill, the coherent scattering depends on the isoscalar scalar, fg, and 

vector, fy, parameters. The latter is effectively multiplied by the average velocity ^ of the 
LSP due to its Majorana nature. These parameters depend on the SUSY model considered. They 
were evaluated in the allowed SUSY parameter space of Kane et ah 0. The construction of the 
dominant scalar parameters suffers from additional uncertainties, which involve the step of going 
from the quark to the nucleon level. In other words, the results are very sensitive to the presence 
of quarks other than u and d in the nucleon. Three such choices indicated by A, B, C are presented 
in Table I. 

From a nuclear structure input point of view, the coherent scattering does not depend on the 
details of the nuclear wave function. It does, however, depend on the nuclear density, i.e. the 
assumed form factor, for fairly massive LSP and correspondingly heavy nuclei. This is because 
the momentum transfer in such cases can be quite high (by nuclear standards) even though the 
energy transfer is small. The form factors used were as realistic as possible [^. The inclusion 
of the form factor results in sizable retardation of the cross section which, for case can be 
a factor of 3 for Ge, and a factor of 14 for Pb. This conclusion is in agreement with previous 
estimates although they do not include detailed calculations for heavy nuclear targets 


32 


B. Spin Contribution 


As we have mentioned in the introduction, the relative importance of the spin matrix element 
is expected to decrease as one goes to the heavier systems vis a vis the coherent scattering. We 
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have seen, however, that the increase due to the mass number in the coherent scattering is offset 
by the decline due to form factor. Thus the spin contribution may not be completely negligible. 

For the spin matrix contribution of we find that its ground state (a | state), in a Ih 

and 2h-lp model space, is more than 95% a 2pi/2 neutron configuration. We have evaluated the 
spin matrix elements up to terms linear in the small amplitudes. Out of the 250 components, only 
two of them contribute to the A = 0 multipole while two more contribute to the A = 2 multipole. 
Thus, we find that, the isoscalar static matrix element suffers from very little retardation, while 
the isovector matrix element is reduced by 17%. Since, in the SUSY parameter space we have 
used, the isovector coupling f\ is larger, this results in a similar retardation of the total matrix 
element. The isoscalar matrix element has other uncertainties. In passing from the quark to the 
nucleon level the amplitude must be multiplied by a factor which ranges from 1, in the Naive 
Quark Model (NQM), to 0.12, if extracted from the EMC data. In our calculations we considered 
both of these extremes. Also, since the dominant configuration is of the neutron variety, the 
isoscalar tends to subtract from the isovector, but this is not so dramatic in our case since the 
isoscalar is smaller in absolute value, especially in the EMC case. 

Our spin matrix elements at g = 0 are 


Oi(0) 


= 0.231 and 


Oo(0) 


= 0.301 


These are a bit larger than the values 


n 2 


Oi(0) ^ Oo(0) 


= 0.2 


29 


Si) 


extracted from Fig. 3 of Ref. [^. They are, however, smaller than the values 


Oi(0) 


n 2 


= 1.00 


and 


Oo(0) 


n 2 


= 1.03 


,9- 

'2 


73 


Ge) 


extracted from Fig. 3 of Ref. . 


As it has been found in Eq. (p!3D, the momentum dependence of the spin part can be described 
in terms of three integrals /qo, loi, hi, where the subscripts indicate the isospin channels. In our 
case, these integrals receive contributions from two multipoles, A = 0 (spin monopole) and A = 2 
(spin-quadrupole). They were judiciously normalized to unity at g = 0. When so normalized these 
integrals are approximately equal (see Fig. 3 and also Ref. |^). We notice that, the monopole 
contribution falls with momentum transfer, as expected, quite fast in fact for large LSP mass. On 
the contrary, the quadrupole contribution starts out at zero and keeps increasing up to about 90 
GeV. As a result its contribution in this mass regime is crucial, since it tends to partly compensate 
for the suppression of the monopole term. We expect these trends to persist even in the most 
elaborate calculations which include the Giant Quadrupole Resonance (GQR) and are currently 
under way. 

As a consequence the momentum transfer suppression is small (a factor of about 4 for mi = 
lOOGeU). Thus, the above two matrix elements take the values 0.120 and 0.153 for the purely 
isovector and isoscalar channels. From Table IV of Ref. one extracts at ?/=0.174 the values 
0.335 and 0.402. Thus the spin matrix element in 207pb 

is less than three times smaller 
compared to that of Ref. for On the other hand the corresponding matrix elements for 
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the A = 29 and A = 207 are similar. When it comes to the cross section the spin contribntion in 
the case of the A = 207 and A = 73 may be similar dne to the presence of rj in Eq. (0), which 
favors the A = 207 system by a factor of abont 2.6 for mi = lOOGeV. 


C. Convoluted rates 


For the coherent cross section in Table III we present resnlts for the parameters and h for 

three typical and experimentally interesting nnclei, Ca, Ge and Pb. From this table we see that, 
the event rate is highest when the intermediate particles and the LSP are lightest (mi = 27GeV, 
case i^2). We notice that, even within the allowed parameter space, the event rate may vary by 
two orders of magnitnde. We also notice that, with the possible exception of the not so realistic 
model A, the Higgs contribntion becomes dominant. This is even more so in models where qnarks 
other than n and d are present in the nnclens with appreciable probabilities, dne to their large 
masses. In the most favorable cases, one may have more than 140 events per year per kilogram 
of Pb target. We notice that the amplitnde dne to the Earth’s annnal motion depends on mi and 
A. In Pb it is less than ±5%. 


Finally, we shonld mention that, for cases and #3 (massive LSP), the rate dne to the 
nnclear form factor can be rednced by a factor of approximately 6. Becanse the form factor 
dependence is more prononnced in heavy nnclei, if the LSP tnrns ont to be massive, the beneht 
of going to heavier targets is somewhat diminished. 


The detection rates obtained, after folding the spin indnced cross section, for are shown 

in Table IV. The folding makes the spin dependence even less sensitive on the momentnm transfer. 
The event rate, however, remains mnch smaller than the coherent scattering for models B and C 
and even for model A in the case ^(^2. These resnlts may snggest that, in the SUSY parameter 
space given in Ref. |^, the spin indnced LSP-nnclens scattering may not be detectable. The 
spin contribntion, however, maybe somewhat enhanced in models in which the LSP is primarily 
a Higgsino . 


In order to estimate the total effect of the nnclear density on the event rates, in Tables III and 
IV the resnlts by ignoring the nnclear form factor (WNFF) are shown. 


VI. SUMMARY AND CONCLUSIONS 


In the present work we have performed calcnlations of the event rate for LSP-nnclens scattering 
for three typical experimentally interesting nnclei, i.e. Ca, Ge and Pb. The three basic ingredients 
of onr calcnlation were the inpnt SUSY parameters, a qnark model for the nncleon and the 
strnctnre of the nnclei involved. 

The inpnt SUSY parameters were calculated in a phenomenologically allowed parameter space 
(cases i^2, #3 of Table I) as explained in the text. In going from the quark to the nucleon 
level, the quark structure of the nucleon was essential, in particular its content in quarks other 
than u and d. For the scalar interaction we considered three models (labeled A, B, C in Table 
I) as described in the text. For the isovector axial coupling one encounters the so-called nucleon 
spin crisis. Again we considered two possibilities depending on the assumed portion of the nucleon 
spin which is due to the quarks (indicated by EMC and NQM in Table I) as described in the text. 

As regards nuclear structure, we employed as detailed as feasible nuclear wave functions. For 
the coherent part (scalar and vector) we used realistic nuclear form factors. For the system 
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we also computed the spin matrix element. The ground state wave function was obtained by 
diagonalizing the nuclear Hamiltonian in a 2h-lp space which is standard for this doubly magic 
nucleus. The momentum dependence of the matrix elements was taken into account and all 
relevant multipoles were retained (in this system one encounters only two, A = 0 and A = 2, due 
to selection rules). 

From our discussion in the previous section we conclude that, even though the spin matrix 
elements are a factor of 3 smaller than those for obtained in Ref. , their contribution 
to the cross section is almost the same (at least for LSP masses around lOOGeV). 

Finally, the obtained results were convoluted with a suitable Maxwell-Boltzmann velocity 
distribution of the LSP’s. This convolution was necessary to partially neutralize the form factor 
retardation, but mainly to compute the modulation of the event rate due to the Earth’s motion. 
We hnd that in almost all cases the event rate due to the Earth’s revolution around the sun is 
very small, less than ±5% around its average value. The event rates thus obtained (see Tables III 
and IV) are, unfortunately, sensitive to the input parameters. 

The inclusion of the nuclear form factor signihcantly retards the event rates for heavy nuclei 
[A > 100) and fairly massive LSP (mi > lOOGeV). However, this retardation does not outweigh 
completely the advantages of using a heavy target. For the spin matrix elements the form factor 
retardation of the usual A = 0 multipole is partially neutralized by the enhancement of the A = 2 
multipole. 

From the data of Tables HI and IV we see that it is possible to have detectable rates ( > 20 
per kilogram per year) only for case ^2 and the realistic nucleon models B and C, associated with 
the scalar Higgs exchange term. In all other cases, including the spin contribution, the calculated 
event rates are too small. 
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APPENDIX 


For the spin contribution of we have to calculate the matrix elements of the operator 


r" = V4^jA(gr)[v\q) 


a 


(72) 


where k = 1 and A = 0, 2 for the states (see Eqs. (^), (^)) 


|(2pi/2) >, 2pi/2) J;j2;l/2> 


(73) 


or the matrix elements 


ME 


'-'0 


± <(2pv2)-1|Tl|(2pv2)-^> + 


(74) 
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± 2^;: Qy, < (2pv2)-‘||r-||(2p./2)-‘[jr‘i2]«:;l/2> 

hh 


The + sign is for isoscalar and the 


< ji||r^||j 2 > are given in ref. 


for isovector matrix elements. The reduced matrix elements 
The relevant radial matrix elements < 'ni/i|j;(gr)|n 2/2 > for 


harmonic oscillator basis can be witten in the compact way 


< nili\ii{qr)\n2l2 > 


^max 

K =0 


X = {lif/i 


(75) 


where 


Kmax = ni + n 2 + m, m = (/i + /2 - 0/2, 


\h ~ 0 | ^ ^ ^ 0 + 0 


The coefficients n2/2, 0 are given by 


^K. 


7rni!n2! 

-4r(?7-i + Zi + |)r(n2 + 0 + §)■ 


ni 712 

E E-!A ^^{nili)K^^{n2l2)K{nl) 


(76) 


where 


n = K,i + K 2 + m, 


A^{nl) 



/'n + / + l/2\ 

V n — k ) 


(see also Ref. |^) and 


0 = 


0 , 

K — ra — 712, 


K — 777 — 772 < 0 
K — 777 — 772 > 0 ’ 


0 , 

K — m — Ki, 


K, — m — Ki < 0 
K — m — Ki > 0 


For the coefficients used in the present work see Ref. |p5 
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Figure Captions 


Fig. 1. 


The integrals Jq (see Eq. (0))> which describe the dominant scalar contribntion (coherent 
part) of the total cross section as a fnnction of the LSP mass (mi), for three typical nuclei: Ca, 
Ge and Pb. The value a/< > = 10“^ was used. 


Fig. 2. 


(a) Plot of the integrals In as a function of the LSP mass mi. This integral gives 
contribution to the LSP-nucleus total cross section for Pb (for the dehnition see Eq. 
the model described in sect. III. The integrals Iqq and /qi are similar, (b) Plot of the 
Iii{u) and Io{u) for Pb, where u is given by Eq. (p!8[) . Note that In is quite a bit less 
compared to Jq. 


the spin 

(H) in 

integrals 

retarded 


Fig. 3. 


Contributions of K integrals (for 1=0 and 1=1) dehned in Eqs. 
rate due to earth’s revolution around the sun: K! 


-(M) and entering the event 


3(c). The other integrals K^ 


and iLgi 


Q in Fig. 3(a), Kq in Fig. 3(b) and K[^ in Fig. 


are similar to K 
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TABLE I. The parameters Pfv^fsyfA’fX fv/ fvy fs ^t'r three SUSY solutions (see Ref. 
113 ). The value of /? = 10 ^ was used. For the definition of fg and fg (models A, B and C) see 
Ref. [jl^ and for the values of tanP, rrih, mn and tua employed see Ref. 0. 


Quantity 

Solution 

Solution #2 

Solution 41^?, 

f3fUZ) 

0.475 X 10-5 

1.916 X 10-5 

0.966 X 10-5 


-1.153 

-1.153 

-1.153 

Pf^q) 

1.271 X 10-5 

0.798 X 10-5 

1.898 X 10-5 

ffiq)/fviq) 

0.222 

2.727 

0.217 

fsiq)/PfUq)i^odelA) 

6.3 X 10-3 

3.6 X 10-3 

2.4 X 10-3 

fsiq)/PfUq)i^odelB) 

0.140 

3.5 X 10-2 

5.8 X 10-2 


1.746 X 10-5 

2.617 X 10-5 

2.864 X 10-5 

fv/fv 

-0.153 

-0.113 

-0.251 

LSP mass{GeV) 

126.0 

27.0 

102.0 

tanl3 

10.0 

1.5 

5.0 

tan2a 

0.245 

6.265 

0.528 

mh,mH,mA 

116,346,345 

110,327,305 

113,326,324 

fg{H) (model A) 

1.31 X 10-5 

1.30 X 10-^ 

1.38 X 10-5 

ff/fs (model A) 

-0.275 

-0.107 

-0.246 

fg{H) (model B) 

5.29 X 10-^ 

7.84 X 10-3 

6.28 X 10-^ 

fliH) (model C) 

7.57 X 10-^ 

7.44 X 10-3 

7.94 X 10-^ 

mz) 

- 

- 

- 

f\{Z) 

1.27 X 10-2 

5.17 X 10-2 

2.58 X 10-2 

flmNQM) 

0.510 X 10-2 

3.55 X 10-2 

0.704 X 10-2 

fAmEMC) 

0.612 X 10-3 

0.426 X 10-2 

0.844 X 10-3 

f\i.q) 

0.277 X 10-2 

0.144 X 10-2 

0.423 X 10-2 

ll(NQM) 

0.510 X 10-2 

3.55 X 10-2 

0.704 X 10-2 

ll(EMC) 

0.612 X 10-3 

0.426 X 10-2 

0.844 X 10-3 


1.55 X 10-2 

5.31 X 10-2 

3.00 X 10-2 


TABLE II. The quantity go (forward momentum transfer) in units of fm ^ for three values of 
nil and three typical nuclei. In determining go the value yj< > = 10“^ was employed (see Eq. 

El))- 


<?o {fm Q 

Nucleus 

mi = 30.0 GeV 

mi = 100.0 GeV 

mi = 150.0 GeV 

Ca 

0.174 

0.290 

0.321 

Ge 

0.215 

0.425 

0.494 

Pb 

0.267 

0.685 

0.885 
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TABLE III. The quantity y ^Kg ^ and the parameter h (oscillation due to the earth’s 

motion around the sun) for the coherent vector and scalar contributions. For the dehnition of 
A, B, C, see text. NFF and WNFF stand for Nuclear Form Factor and Without Nuclear Form 
Factor, respectively. 




Vector Contribution 

Scalar Contribution 





h 

/ dN 
\ dt 


) 

h 


Solution 



Model A 

Model B 

Model C 



#1(AFF) 

0.264 X 10-3 

0.029 

0.151 X 10-3 

0.220 

0.450 

-0.002 


#1(IFAFF) 

0.378 X 10-2 

0.102 

0.212 X 10-2 

3.087 

6.322 

0.054 


#2 (AFP) 

0.162 X 10-3 

0.039 

0.410 X 10-1 

142.860 

128.660 

0.026 

Pb 

#2(IFAFF) 

0.332 X 10-3 

0.102 

0.872 X 10-1 

303.390 

273.220 

0.054 


#3 (AFP) 

0.895 X 10-3 

0.038 

0.200 X 10-3 

0.377 

0.602 

-0.001 


#3(IFAFF) 

0.970 X 10-2 

0.102 

0.220 X 10-2 

4.114 

6.576 

0.054 


#1(AFF) 

0.151 X 10-3 

0.043 

0.779 X 10-1 

0.120 

0.245 

0.017 


#1(IFAFF) 

0.518 X 10-3 

0.102 

0.230 X 10-3 

0.353 

0.723 

0.054 

Ge 

#2(AFF) 

0.527 X 10-^ 

0.057 

0.146 X 10-1 

51.724 

46.580 

0.041 


#2(IFAFF) 

0.708 X 10-^ 

0.102 

0.196 X 10-1 

69.506 

62.595 

0.054 


#3 (AFP) 

0.481 X 10-3 

0.045 

0.101 X 10-3 

0.198 

0.316 

0.020 


#3(IFAFF) 

0.144 X 10-2 

0.102 

0.267 X 10-3 

0.525 

0.839 

0.054 


#1(AFF) 

0.790 X 10-^ 

0.053 

0.340 X 10-1 

0.055 

0.114 

0.037 


#1(IFAFF) 

0.135 X 10-3 

0.102 

0.520 X 10-1 

0.085 

0.174 

0.054 

Ca 

#2 (AFP) 

0.264 X 10-3 

0.060 

0.612 X 10-2 

22.271 

20.056 

0.048 


#2(IFAFF) 

0.307 X 10-3 

0.102 

0.704 X 10-2 

25.601 

23.055 

0.054 


#3 (AFP) 

0.241 X 10-3 

0.053 

0.435 X 10-1 

0.090 

0.144 

0.038 


#3(IFAFF) 

0.388 X 10-3 

0.102 

0.642 X 10-1 

0.133 

0.212 

0.054 


TABLE IV. The spin contribution in the LSP-nucleus scattering in for two cases: EMC 

data and NQM Model (see sect. II). For other notations see caption of Table III. 



EMC DATA 


NQM MODEL 

Solution 

(f )„ iv-'Kg-') 

h 

(f )„ iv-'Kr') 

h 

#1(AFF) 

0.285 X 10-2 

0.014 

0.137 X 10-2 

0.015 

#1(IFAFF) 

0.130 X 10-1 

0.054 

0.551 X 10-2 

0.054 

#2(AFF) 

0.041 

0.046 

0.384 X 10-2 

0.056 

#2(IFAFF) 

0.062 

0.054 

0.405 X 10-2 

0.054 

#3(AFF) 

0.012 

0.016 

0.764 X 10-2 

0.017 

#3(IFAFF) 

0.047 

0.054 

0.269 X 10-1 

0.054 
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